Substructure lensing in galaxy clusters as a constraint on low-mass sterile neutrinos 
in tensor- vector-scalar theory: The straight arc of Abell 2390 
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Certain covariant theories of the modified Newtonian dynamics paradigm seem to require an additional hot 
dark matter (HDM) component - in the form of either heavy ordinary neutrinos or more recently light sterile 
neutrinos (SNs) with a mass around 1 leV - to be relieved of problems ranging from cosmological scales down 
to intermediate ones relevant for galaxy clusters. Here we suggest using gravitational lensing by galaxy clusters 
to test such a marriage of neutrino HDM and modified gravity, adopting the framework of tensor- vector-scalar 
theory (TeVeS). Unlike conventional cold dark matter (CDM), such HDM is subject to strong phase-space con- 
straints, which allows one to check cluster lens models inferred within the modified framework for consistency. 
Since the considered HDM particles cannot collapse into arbitrarily dense clumps and only form structures well 
above the galactic scale, systems which indicate the need for dark substructure are of particular interest. As 
a first example, we study the cluster lens Abell 2390 and its impressive straight arc with the help of numer- 
ical simulations. Based on our results, we outline a general and systematic approach to model cluster lenses 
in TeVeS which significantly reduces the calculation complexity. We further consider a simple bimodal lens 
configuration, capable of producing the straight arc, to demonstrate our approach. We find that such a model is 
marginally consistent with the hypothesis of 1 leV SNs. Future work including more detailed and realistic lens 
models may further constrain the necessary SN distribution and help to conclusively assess this point. Cluster 
lenses could therefore provide an interesting discriminator between CDM and such modified gravity scenarios 
supplemented by SNs or other choices of HDM. 

PACS numbers: 04.50.Kd, 98.80.-k, 02.60.Lj 



I. INTRODUCTION 

The current concordance model of cosmology has been 
remarkably successful in forming a consistent picture on 
the largest physical scales. For instance, it provides suit- 
able explanations for the observations on supernovae la fl], 
large-scale structure J2] |3), weak lensing J4), and the cos- 
mic microwave background (CMB) J5] |6). Based on gen- 
eral relativity (GR), this so-called ACDM model relies on two 
phenomenologically motivated constituents, cold dark mat- 
ter (CDM) and dark energy (DE), which account for approx- 
imately 95% of the energy-matter content of the Universe. 
From a more fundamental point of view, however, the intro- 
duction of a DE component is theoretically challenging and 
extremely fine-tuned, despite the many proposals for its dy- 
namics [7|. On the other hand, the concept of CDM also suf- 
fers from several issues such as the lack of direct experimental 
detection [8 |, the question of its cosmological abundance ||9l , 
and problems related to the formation of structure on small 
scales CM!- 

A perhaps more natural solution might be that gravity gen- 
uinely differs from GR, which expresses itself as either one 
or even both of the above dark components |[T4HT9l . Here 
we want to consider a particular class of modifications which 
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give rise to what has become known as the modified Newto- 
nian dynamics (MOND) paradigm [20, 21 1. On a nonrelativis- 
tic level, MOND aims at solving the missing mass problem 
by postulating an acceleration-dependent change of Newton's 
law which is characterized by a scale «o l22l : 

pi — )a = -V® N + S. (1) 
\ao) 

Here, denotes the common Newtonian potential of a mat- 
ter source and S is a solenoidal vector field determined by the 
condition that a can be expressed as the gradient of a scalar po- 
tential. The function fi, controlling the modification of New- 
ton's law, has the following asymptotic behavior: 

u(x) ~ x x <s 1 , 

- 1 , (2) 

H(x) ~ 1 x » 1 . 

Equation ([T| has been constructed to agree with the fact that 
the rotation curves of spiral galaxies become flat outside their 
central parts. Analyzing observational data, Milgrom esti- 
mated an acceleration scale of ao » 1.2 X 10~ 10 m s~ 2 . 

The MOND paradigm still appears suitable to explain the 
observed "conspiracy" between the distribution of baryons 
and the gravitational field in spiral galaxies E3l f2M . It is 
striking that such a simple prescription leads to extremely suc- 
cessful predictions for galaxies ranging over five decades in 
mass (see, e.g., Refs. [25 , 26] for reviews), including our own 
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Milky Way ||27ti29l . dwarf spheroidals B0ij32l . x-ray dim el- 
liptical galaxies 11331 [3411 . and tidal dwarf galaxies ll35l [361 . 
In addition, MOND successfully reproduces empirical galaxy 
scaling relations such as the well-known Tully-Fisher relation 
07113811 . and more recently the central surface brightness pre- 
dicted by dark halos I139WT1 . 

While the framework of MONDian dynamics works ex- 
tremely well on galactic scales, the situation in galaxy groups 
and clusters is quite different: Several studies of such systems 
Il42lf44l have shown that an additional nonluminous matter 
component is required to explain observations, even after tak- 
ing into account the gravitational boost induced by the MOND 
formula. Assuming that MOND is a viable description for 
such gravitating systems, this result has led to the question 
of what the needed matter component should be. It is obvi- 
ous that any possible form of exotic CDM is disfavored as it 
would cause the original idea of Eq. ([TJ to become redundant. 

Another problem arises from the fact that the original 
MOND formulation does not specify cosmology or the nature 
of gravitational light deflection. Recent developments in the 
theory of gravity, however, have been able to embed MON- 
Dian dynamics into fully Lorentz-covariant theories by means 
of a dynamical four- vector field B3H491 . Although still lack- 
ing a derivation from fundamental principles underpinning the 
MOND paradigm, these theories allow for new predictions re- 
garding cosmology and structure formation 1 50-53] as well 
as gravitational lensing [54-60]. An appealing feature of such 
modifications is also that they might be able to simultaneously 
explain the observed effects of DE [18 61-66 1, but we do not 
consider this possibility in the present work. 

Adopting the framework of tensor-vector-scalar theory 
(TeVeS) B31l . it was possible to investigate the resulting CMB 
and matter power spectra 11501 . It has been found that the 
known baryonic matter content of the universe is not enough 
to reconcile the theory with observations fl67l . To overcome 
this problem, one may resort to a model with ordinary neutri- 
nos of mass around 2eV which is able to describe the obser- 
vational data in a qualitatively acceptable way, but the cor- 
responding fits do not match the excellent agreement of a 
ACDM model. Interestingly, the idea of massive neutrinos 
has already been discussed to provide a solution to the lack of 
matter on cluster scales Il42l[43l , and to explain the observed 
weak lensing map of the galaxy merger 1E0657 - 558 ("bul- 
let cluster") [57 , 68 1. It should be mentioned that the needed 
neutrino mass of 2eV is barely consistent with the current up- 
per limit on the electron neutrino's mass (e.g., see Ref. [69 1); 
future measurements such as the Karlsruhe Tritium Neutrino 
Mass Experiment (KATRIN) 0(3 ED will be able to explore 
a mass range well below the 2eV threshold. 

Nevertheless, the rather unsatisfactory results of this so- 
lution on large scales, especially for the CMB anisotropy 
power spectrum, and problems within galaxy groups tl44ll have 
led to deem the hypothesis of very massive ordinary neutri- 
nos unattractive. Alternatively, the required additional matter 
could be provided in the form of heavy (right-handed) sterile 
neutrinos (SNs) which are motivated by theoretical consider- 
ations in particle physics (e.g., see Refs. ir72]474l and refer- 
ences therein) and offer an elegant way to explain the small 



masses of active neutrinos via the "seesaw mechanism" 11751 - 
l77l . The conceptual advantage of such an approach lies in 
combining the success of modified gravity on small scales 
with new physics in a sector of the standard model which 
is known to be incomplete ll78l and in need of revision IT791 . 
Motivated by a possible interpretation of the MiniBooNE ex- 
periment (80), Angus [81 1 has suggested to use a single light 
species of SNs with a mass of approximately 1 leV and inves- 
tigated its consequences. If such SNs decouple while they are 
relativistic and in thermal equilibrium, one should obtain both 
a background evolution and a CMB power spectrum which 
are basically indistinguishable from a standard ACDM cos- 
mology [82 1, while at the same time, this additional hot dark 
matter (HDM) component may give rise to a correct predic- 
tion of the linear matter power spectrum and represents a suit- 
able candidate for the missing mass in galaxy clusters without 
spoiling MONDian dynamics on galactic scales 11831 . As for 
the nonlinear regime of structure formation, the situation is 
still unclear. Because of the more sophisticated mathemati- 
cal structure of the nonlinear TeVeS field equations (or that of 
related theories) as opposed to those of GR, there seems cur- 
rently no way to gain reliable information about the nonlinear 
evolution. This difficulty is somewhat reflected by the fact that 
the resulting field equations in the quasistatic, nonrelativistic 
limit typically remain highly nonlinear. Assuming an ad hoc 
modification of the original MOND formula Eq. ([TJ, however, 
a first simplified attempt into this direction is discussed in Ref. 

em- 
it is noteworthy that TeVeS or TeVeS-like theories in com- 
bination with sufficiently abundant massive neutrinos provide 
the most consistent relativistic MOND framework presented 
in the literature so far 11851 : nevertheless, there are still in- 
numerable aspects which need to be tested further. A pow- 
erful astronomical tool to challenge these theories is gravi- 
tational lensing. Using spherically symmetric lens models, 
TeVeS has been tested against the CfA-Arizona Space Tele- 
scope Lens Survey data of image-splitting lens galaxies lf56l . 
building on earlier work in the context of MOND lf55ll (also 
see Refs. [86, 87 1). Later, the analysis has been extended to 
nonspherical lens geometries [88 1. While there is currently 
disagreement on how well TeVeS explains the lensing proper- 
ties of relatively isolated galaxies [86 , 87 89 1, the models con- 
sistently seem to fail for lens galaxies embedded into dense 
environments like groups or clusters. A detailed individual 
discussion on such lens systems can be found in Ref. ||88l . 
Several authors [56, 88 1 have argued that this issue might be 
attributed to the already known dynamical problem in groups 
and clusters and therefore resolved if the actual environment 
of lens galaxies is properly taken into account. Another way 
of testing the theory is offered by weak galaxy-galaxy lens- 
ing. Using data from the Red-Sequence Cluster Survey and 
the Sloan Digital Sky Survey (SDSS), it has been found that 
the most luminous galaxies (> 10 n L o ) would require a sub- 
stantial fraction of nonbaryonic matter [90|. Although this 
result needs to be confirmed by larger data sets before a firm 
conclusion can be drawn, it might hint towards a problem with 
the original MOND idea on galactic scales. Again, SNs with 
a mass around 1 leV could provide a remedy as they should be 
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able to cluster densely enough in such massive systems 
However, it remains to be seen whether such an approach can 
explain observations. 

In this work, we suggest to test TeVeS and the massive SN 
hypothesis in the context of complex lens systems which are 
typically present in the central regions of galaxy clusters. A 
previous analysis [91 1 already revealed that such an environ- 
ment can put stringent constraints on the distribution and plau- 
sibility of the needed dark neutrino component, thus provid- 
ing an excellent testbed for our purposes. Generally, the ad- 
vantage of galaxy clusters lies in the independent estimates of 
baryonic matter, inferred from observed X-ray and stellar lu- 
minosities, and of the system's total mass distribution based 
on a combination of weak and strong gravitational lensing. 
Being insensitive to the dynamical state of the deflecting mass, 
the latter techniques are particularly suited to constrain the 
properties of the dark component. In contrast to weak lensing 
estimates, strong lensing is basically free of statistical uncer- 
tainties and offers a unique and robust probe of the matter dis- 
tribution on scales < lOOkpc. Here we shall use strong lensing 
to further test the viability of 1 leV SNs. Unlike conventional 
CDM, light SNs are subject to strong phase-space bounds set 
by the Tremaine-Gunn limit ||92ll , which allows one to check 
cluster lens models inferred within the modified framework 
for consistency. Since this limit prevents SNs from clustering 
into dense clumps, galaxy cluster lenses with a considerable 
amount of dark substructure provide an ideal target for our in- 
tentions. As a first example, we shall study the galaxy cluster 
Abell 2390 (A2390) with its notorious straight arc, and inves- 
tigate whether it is possible to reproduce this particular lens 
feature in TeVeS. 

The paper is structured as follows: Starting with a brief in- 
troduction to TeVeS and the formalism of gravitational lens- 
ing in Sec. [II] we give an observational summary of the galaxy 
cluster A2390 and its pronounced straight arc in Sec. [Ill] Con- 
tinuing with the setup for a simplified density model of A2390 
in Sec. |IV| we discuss results for quasiequilibrium configura- 
tions in Sec. [V] Based on the latter, we outline a systematic 
approach to cluster lenses in TeVeS, and describe a lens model 
for the straight arc in Sec. VI Finally, we conclude in Sec. 



VII For clarity, technical and numerical details are given in 
an appendix. 



paradigm. Compared to GR, TeVeS is based on the well- 
known Einstein metric g^ v and two additional dynamical 
(gravitational) fields, a four-vector and a scalar (/>; it also 
invokes a second "physical" metric g^ v , related to g^ v through 
a nonconformal transformation involving the new fields, and 
being exclusively used for the coupling of matter fields (p, v — 
0,...,3): 



-2U„U> 



sinh(20). 



(3) 



Furthermore, the theory introduces a new set of free parame- 
ters: The constants k and K as the couplings of the scalar and 
vector field, respectively, and a length scale I. Albeit not a 
unique extension, TeVeS is the most popular "MONDian rep- 
resentative" so far, and a variety of its aspects have been ex- 
tensively studied in the literature (see Ref. 11931 for a review). 
Although the original formulation of TeVeS suffers from sev- 
eral problems, e.g. in the strong gravity regime [94, 95 1 or - at 
least for certain models - in the cosmological domain 15311961 , 
it still provides a viable description of relativistic MOND on 
extragalactic scales. 

In what follows, we shall restrict ourselves to weak fields 
and quasistatic systems l97l . Under these premises, the usual 
metric potential <t> approximately takes the form 



<J> SB + (p, 

where the TeVeS scalar field (p obeys 

V • [n (kl 2 (V(p) 2 ) Vcf>] = kGp, 



(4) 



(5) 



and 0# corresponds to the common Newtonian potential 
given by Poisson's equation. In order to reasonably recover 
the dynamics of MOND, the constants k and K have to be 
much less than unity, i.e. k,K<zz 1, and G ~ Gn, but one 
should also see Ref. RBI for a discussion on lower bounds 
of k. The function fi(y), which is related to a potential-like 
term in the scalar field's action, has to satisfy the following 
conditions: 



j"O0 
H(y)~ 



h y 

y 



b' 



y « l, 



(6) 



II. TEVES AND GRAVITATIONAL LENSING 
IN A NUTSHELL 



In the following, we will briefly review TeVeS and gravita- 
tional lensing within the approximations used for quasistatic 
systems like galaxies and galaxy clusters. Henceforth, we 
shall use units with c = 1, and the "square" of a vector de- 
notes the square of its Euclidean norm, i.e. A 2 = ||A||2 2 . 



A. The nonrelativistic limit of TeVeS 

TeVeS theory [45 ] is a fully relativistic description of grav- 
ity which provides a covariant framework for the MOND 



where 



y = kl 2 (Vcf>) 2 , 



(7) 



and b is a positive real constant. The TeVeS length scale / can 
be related to Milgrom's ao by 



a 



4nZl ~ Anl ' 



(8) 



where H = 1 - K/2 - 2(p c , and <p c is the cosmological value of 
the scalar field which is assumed to be small (<p c <sz 1) |45l . If 
not specified in any other way, we shall work with 



b = l, k = 0.01 



(9) 
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throughout this paper. This is justified following the analysis 
of Ref. [59 1 : While the value of b is absorbed into the defini- 
tion of the scale «o in Eq. ([8j, the TeVeS lensing maps have 
been shown to be generally insensitive to variations of the pa- 
rameter k as long as it is small, k < 0.01. A particular choice 
for the free function /j(y), suitable for the gravitational lensing 



analysis of A2390, will be given and discussed in Sec. IV 



B. Background cosmology 

Assuming the basic fields to partake of the symmetries 
of the Friedmann-Robertson-Walker (FRW) spacetime, the 
TeVeS analog of Friedmann's equation in the Einstein frame, 
delineated by the metric g^, reads 



8ttG 



A 
3' 



(10) 



where we have included a cosmological constant A, and i; 
the energy density of the scalar field. Using Eq. ([3]l for g^ v 
we switch to the matter frame and finally obtain 



with 



1 da 
a dt 



dt 



'(;-*)• 



(11) 



(12) 



According to previous studies |45, 50], it is consistent to as- 
sume that the cosmological scalar field evolves slowly in time 
throughout cosmological history. Thus, its contribution to the 
Hubble expansion is negligibly small, with a ratio O(k) com- 
pared to the matter contribution. Setting = and recalling 
that <p <k 1 at the background level, the physical Hubble pa- 
rameter can be expressed as 



H 2 x H 2 x Hi (Q„,(l + z) 3 + O a + £i*(l + z) 2 ) , 



(13) 



where = 1 — Q A - Q m for a flat universe. 

To calculate angular diameter distances in the context of 
gravitational lensing, we further need to specify a particular 
choice of parameters. Based on the assumption of a single 
species of lleV SNs llSTl 1851 . we shall use a flat cosmology 
with Q. m = Q. h + Q v = 0.29 and h = 0.7, where the Hub- 
ble constant is Ho = 100/;km s _1 Mpc -1 . Note that this gives 
a background evolution which is virtually indistinguishable 
from a standard ACDM model. 



C. Light bending in TeVeS and lensing formalism 

Being a metric theory, the propagation of light rays in 
TeVeS is determined by the null geodesies of the physical 
metric g^. If the metric potential <$> given by Eq. Q and 
the peculiar velocity v of the lens are small (<I>, v <s 1), one 
can presume a locally flat spacetime being disturbed by the 
potential O; these conditions are well satisfied for galaxies 



and galaxy clusters. Following the lines of Ref. [98 1, the de- 
flection angle of a light ray under these assumptions can be 
expressed as 



a 



x Odl = a G R + 2 I V±$dl. 



f 



(14) 



where Vj. denotes the two-dimensional gradient operator per- 
pendicular to light propagation, and integration is performed 
along the unperturbed light path (Born's approximation). In 
addition to the deflection angle caused by the Newtonian po- 
tential (D;v, there is a contribution arising from the scalar field 
<p. Assuming that the metric potential O is known from solv- 
ing the corresponding field equations, one can directly pro- 
ceed to calculate the usual lensing quantities, fully adopting 
the standard GR formalism P31I591 which is briefly reviewed 
below. 

In gravitational lensing, it is convenient to introduce the de- 
flection potential Y(0) |§3S): 

D ds 



W(d) = 2- 



I 



®(D d 0,z)dz, 



(15) 



D s D d 

where we have used 9 = $,/D d and have chosen coordinates 
such that unperturbed light rays propagate parallel to the z 
axis. Here £ is the two-dimensional position vector in the 
lens plane, and D s , D d , and D ds are the (angular diameter) 
distances between source and observer, lens and observer, and 
lens and source, respectively. If a source is much smaller than 
the angular scale on which the lens properties change, the lens 
mapping can locally be linearized. Thus, the distortion of an 
image can be described by the Jacobian matrix 



dJ3_ 
36 



K-y\ 

-72 



1 



-72 

k + 7\ 



(16) 



where (3=r]/D s and r\ denotes the two-dimensional position 
of the source. The convergence k is directly related to the 
deflection potential through 



1 



1 ( d 2x ¥ d 2x ¥ 
d6 2 dOj 



2 " 2 

and the shear components j\ and j2 are given by 
\(d 2x Y d 2x V\ d 2y ¥ 



(17) 



7i 



del 



de\ j 



72 



<96>i<9<9 2 



r : 



Points in the lens plane where 



det^l = (1 - k- y)(l - k + y) = 0, 



(18) 



(19) 



form closed curves, the critical curves. Their corresponding 
image curves residing in the source plane are called caustics. 
Images near critical curves can be significantly magnified and 
distorted, which, for instance, is indicated by the giant lumi- 
nous arcs formed from source galaxies near caustics. Knowl- 
edge about the exact shape and location of these curves al- 
ready allows one to make solid statements about the system's 
strong lensing properties. Since the possible critical curves 
and caustics of A2390 are relatively well constrained by ob- 
servations of the here considered straight arc, estimating them 
for several density models in TeVeS will be one of the main 
tasks in this work. 



5 




FIG. 1. A small section of an HST/WFPC2 observation of A2390 
shows the impressive straight arc on the left side. Characterized by 
two breaks along its light profile (present in other observed bands as 
well), the arc can be decomposed into three segments labeled A, B 
and C, respectively |99|. Also visible are the galaxy 2592, which is 
located adjacent to the arc, as well as the galaxy 6666 (see Table[I). 



III. OBSERVATIONS OF THE GALAXY CLUSTER A2390 
A. X-ray gas and member galaxies of A2390 

The galaxy cluster A2390 at redshift z = 0.23 HUO] [TUTl 
is one of the richest and most luminous clusters known in the 
literature. Several interesting properties, e.g. the large abun- 
dance of lensing arcs and arclets j99l, an elongated galaxy 
distribution 1 102 1 and its large velocity dispersion 1 103), have 
made the analysis of this system particularly attractive. In the 
context of GR, A2390 has been subject to extensive study 
by means of different techniques including virial (e.g., Ref. 
1110410 . x-ray II105H108I . redshift-space caustic 111091 and both 
weak IfTTOlJPHl and strong HI 15tfTT8l lensing studies. 

Observations with CHANDRA exhibit a very concentrated 
and highly peaked x-ray emission, indicating a strong cooling 
flow which is centered on cluster's central cD galaxy 11061 . 
On large scales, the x-ray morphology has been found to be 
strongly elliptical with an overall position angle (PA) compa- 
rable to the main cluster direction in the optical (PA = 133°) 
1991 . Here and below, the PA is defined as the angular offset of 
the major axis with respect to the north-south direction, being 
measured counterclockwise. The data provide evidence for an 
elongated x-ray morphology in the very central part, and sug- 
gest the existence of a substructure in the cluster gas located 
roughly 40" (~ 147kpc) from the cluster center. The CHAN- 
DRA image further reveals large-scale cavities in the x-ray 
surface brightness extending approximately 400kpc from the 
center, where a sharp break in the surface brightness profile is 
visible. As observed in several other clusters 111 191 . such cav- 
ities are likely produced by bubbles of radio plasma emitted 
by the central active galactic nucleus. Despite these irregu- 



larities and the appearance of a secondary gas peak, however, 
the x-ray observations indicate that the system as a whole is 
relatively regular and, to good approximation, dynamically re- 
laxed. Thus, if one excludes the cluster's central part, the over- 
all assumption of hydrostatic equilibrium appears as a reason- 
able one. 

There are also several studies of individual galaxies within 
the cluster. For instance, the properties of the central cD 
galaxy have been examined using optical 111201 . infrared and 
radio observations 11211 . A large sample of 216 confirmed 
cluster members based on photometric and spectroscopic in- 
formation is presented in Ref. II 1 221 . More recent observa- 
tions include a selection of 48 early-type member galaxies 
which has been used to investigate their evolutionary status 
11231 . We note that the available observational data will be 
important for building a realistic cluster model in TeVeS. 



B. The straight arc of A2390 

Among several arcs and arclets, the cluster A2390 exhibits 
an unusual, strongly lensed straight arc (see Fig. [TJ which is 
located approximately 38" (~ 140kpc) from the central cD 
galaxy [99 1. This particular arc is unusual in the sense that, 
as it is both located in the outer core region and adjacent to a 
lens galaxy lying in between arc and cD galaxy, it would be 
expected to appear curved with respect to the massive cluster 
center or the closest galaxy. Along its light profile, the arc 
further exhibits two breaks in surface brightness, symmetri- 
cally located relative to the closest galaxy's center. Spectro- 
scopic analysis of the arc revealed that it is actually the joint 
image of two different sources, one at redshift z = 0.913 (cor- 
responding to B - C in Fig. [j} [99] and the other at z = 1.033 
(corresponding to A) II124I . In addition, ISOCAM observa- 
tions of the image segment B - C indicate the presence of an 
active star forming region and support the scenario of two in- 
teracting source galaxies atz = 0.913 11251 . Nevertheless, the 
found straightness requires a rather special lens configuration 
(also see Sec. |V A\ . 

Apart from the system A2390, there also exist other detec- 
tions of (relatively) straight images which are typically well 
modeled from the visible distribution of bright galaxies helped 
by the central cluster potential 1 126 - 128 1. As already pointed 
out in the literature B115I . a similar approach for A2390 within 
the usual framework of CDM would require extremely high 
mass-to-light ratios for individual galaxies, and thus yields a 
rather unrealistic scenario. In recent years, several authors 
have considered possible lens models which aim at reproduc- 
ing such a straight image, and a first attempt was performed 
in Ref. [99|. For instance, the fold caustic of a single, highly 
elliptical cluster lens can be used to create a straight image 
111 161 . Such a model gives a result comparable with the arc's 
morphology, but fails to explain infrared observations. Adopt- 
ing a very large ellipticity of the central cluster profile, it 
was demonstrated how a cusp model may produce the de- 
sired elongated image morphology HI 181 : however, this so- 
lution seems incompatible with other lensing constraints of 
the system. Building on the existence of x-ray substructure in 
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the arc's vicinity, the authors of Ref. B117I employed a two- 
component model using an elliptical cluster center with axis 
ratio b/a = 0.7 to explain the arc. Despite a slight deviation at 
~ lcr significance, the obtained x-ray temperature profile and 
the projected mass within 38" (~ 140kpc) appear consistent 
with those derived from the observed x-ray luminosity II106L 
It seems obvious that any suitable model needs substan- 
tial fine-tuning to form the necessary lens configuration for 
straight images. As a consequence, all of these models are 
extremely sensitive and unstable with respect to perturbations 
due to the closest galaxy or additional substructure in the in- 
tracluster medium (ICM). While this does not pose a problem 
per se, it is nevertheless interesting to look for models with 
improved stability. From a general analysis on how to form 
straight images 11 151 . it has been concluded that the most 
likely configuration involves a dark mirror component of the 
nearest galaxy located on the opposite side of the arc, coun- 
terbalancing the effect of the visible galaxy. With the help of 
the central cluster profile, this yields a so-called beak-to-beak 
model which explains the observed straight arc and, if real- 
ized with such a "dark galaxy," is sufficiently stable against 
local perturbations. Alternatively, there is also the possibility 
of a lips catastrophe 11 151 . i.e. a lips caustic just emerged or 
just about to emerge in three-dimensional caustic space (for a 
demonstration of a lips catastrophe in A370 see Ref. II129I0 . 
Since such a model requires the lensing convergence - equal 
to the projected matter density in GR only - to peak at the arc's 
position, however, it is not supported by observations. 

C. A challenge for TeVeS and hot dark matter 

Concerning the situation in TeVeS, we may already state 
that the "dark galaxy" approach, i.e. a nonluminous mat- 
ter distribution of galactic size, cannot be achieved with our 
choice of lleV SN HDM. Assuming that these particles are 
relativistic and thermalized at the time of decoupling (just 
like for active neutrinos around a temperature of several MeV 
which is much larger than the considered mass) |130|, their 
Fermi-Dirac distribution freezes in, and their phase-space 
density is constrained by the Tremaine-Gunn (TG) limit |92|. 
For instance, a HDM galaxy in TeVeS would have a typical 
phase-space occupation number (we neglect factors of n and 
order unity) 

fv>dN M I h \ 3 a 2 Q h 3 
d 3 xd 3 p m Y \m v crrc ] Gm^cr 5 



VlleV/ \ 100km s" 1 / ' 

which exceeds unity, and thus the TG limit for thermal relics, 
unless the HDM mass m v is much larger than lleV (e.g., 
~ IkeV warm dark matter) and/or the structure's velocity dis- 
persion cr » 100km s , hence above the galactic scale. The 
estimate given in Eq. ( p0] > assumes that the structure's dense 
core is subject to the Newtonian regime (p ~ 1), which gives 
a core size r c ~ GM/cr 2 , and the total mass M ~ cr 4 /Gao 
can be well approximated within the "deep-MOND" limit 



(p ~ -y(y). Also note that moving to masses significantly larger 
than m v = 1 leV would spoil the dynamics of MOND in galax- 
ies and thus eliminate the use of such HDM in the first place 

Therefore, a combination of HDM and modified gravity 
may, in principle, face a challenge in order to create observed 
effects of dark substructure. The TG phase-space bound not 
only applies to HDM substructure, but also to its global dis- 
tribution within the cluster, which presents a well-posed and 
constraining general test of TeVeS or similar theories supple- 
mented by an additional HDM component. As other realistic 
lens models for the straight arc II 1 1 71 also suggest a substan- 
tial amount of dark substructure, a basic question is whether 
there are TeVeS lens models which are compatible with the 
TG bound for lleV SNs. Before we can address this point, 
however, we need a reliable way of modeling the straight arc 
in TeVeS. An approach into this direction will be discussed 
below. 

In preparation for the following sections, we introduce the 
terminology and procedure used for two different kinds of lens 
configurations in our analysis: 

Quasiequilibrium configurations Here we consider con- 
figurations which are based on the assumption of hydrostatic 
equilibrium. Both the cluster gas and the (SN) HDM compo- 
nent are modeled by symmetric, central density distributions, 
the latter having a maximum phase-space density set by the 
TG limit which is inferred in a self-consistent way by con- 
sidering the equation of state for a partially degenerate neu- 
trino gas ||83~1|131I . In addition, we include substructure in the 
form of visible galaxies and further allow for perturbations of 
the central distribution (gas + HDM) which are modeled by 
the same density profile as the central one (corresponding to 
structure of equal scale). We then check whether such config- 
urations can produce the observed straight image in TeVeS. 

Nonequilibrium configurations In this case, we allow for 
any HDM distribution which is capable of explaining the 
straight arc. This includes complex distributions with multi- 
peaked mass densities and concentrations of different scale. 
Although we outline a general approach to lens models in 
TeVeS, we restrict our analysis to a bimodal configuration 
based on a model in GR (cf. Table [II] below) whose com- 
ponents exhibit dispersions cr > 500km s and appear con- 
sistent with the crude estimate of Eq. (GDI, i.e. cr > 
400km s _1 (ot v /1 leV) 4 ' 5 . Approximately treating each den- 
sity peak as a symmetric equilibrium distribution of SNs, we 
investigate whether they satisfy the TG phase-space limit for 
m v - lleV. For simplicity, we do not account for baryonic 
substructure (galaxies) in this context. 



IV. QUASIEQUILIBRIUM MODEL OF A2390 

Because of the nonlinear relation of the TeVeS scalar field 
to the underlying matter distribution, we cannot work with 
projected quantities, but need to perform our calculations in 
three dimensions. This significantly complicates the lensing 
analysis of A2390 and requires knowledge about the clus- 
ter's three-dimensional matter density. A first approach to our 
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TABLE I. Positions, line-of-sight configurations and masses of individual galaxy components for the density model of A2390: At the cluster's 
redshift (z = 0.23), an angular scale of 1" corresponds to approximately 3.7kpc. 
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FIG. 2. TeVeS equilibrium configuration of lleV sterile neutrinos 
in A2390: The figure shows the calculated density distribution of 
neutrinos (dashed line), the analytic fit to this density using the pro- 
file specified in Eq. {27} (solid line) and the central baryonic matter 
distribution derived from x-ray observations (dotted line). 



problem is to consider cluster configurations which are based 
on the assumption of hydrostatic equilibrium. 



A. Distribution of baryonic material 

Using available data of x-ray gas [ 106 108 1 and individual 
galaxies II 1231 . we have modeled the distribution of baryons 
in A2390. Here we shall briefly present the results which are 
relevant for the analysis in Sec. [V] A detailed description of 
our procedure can be found in App. [A] 

Figure [2] shows the density distribution inferred from x-ray 
observations with CHANDRA (dotted line). In addition to 
this central profile, we consider the contribution of five mas- 
sive early-type galaxies which are located close to the straight 
arc. The masses of these galaxies are derived following a 



twofold approach: The first estimate (denoted as M\ ) is based 
on a direct conversion of observed luminosity to stellar mass 
while the second one {Mi) uses a dynamical method. In what 
follows, we shall consider both prescriptions and present re- 
sults for the two different mass estimates below. We further 
assume that all galaxies can be described by a spherical den- 
sity profile which is closely related to the Hernquist profile 
[132] for elliptical galaxies (see App. |A 2) . Using the nota- 
tion of Ref. 111231 , the basic properties of our models for the 
galactic components are illustrated in Table [I] 

B. Adding massive neutrinos 

As mentioned in Sec. [j] TeVeS requires an additional mat- 
ter component to consistently describe observations of galaxy 
clusters. Assuming lleV SNs within the original formula- 
tion of MOND, the authors of Ref. [83] derived their corre- 
sponding equilibrium density and (radial) velocity dispersion 
distributions for a sample of 30 galaxy groups and clusters, in- 
cluding the system A2390. Starting from the observed density 
and temperature of the ICM [108], p x (r) and ksT x (r), respec- 
tively, the assumption of hydrostatic equilibrium immediately 
allows one to determine the gravitational field as a function of 
radius: 

, . -k B T x (r) (d log p x (r) d log k B T x (r)\ 

g(r)= — 7: + j: , (21) 

wm p r \ a log r alogr / 

where w w 0.6 is the mean molecular weight and m p the mass 
of the proton. The such derived result is typically accurate to 
~ 10% if equilibrium is realized 111081 . Using the above, one 
directly obtains the total enclosed MOND mass which is given 
by 

, r 2 g(r)p.(x) g 

M(r) = , x = — . (22) 

G flo 

Here fi corresponds to the MOND interpolating function de- 
fined in Eqs. ([TJ and Note that this is the only stage where 
the modification of gravity is involved. Once this function is 
specified, Eq. |22} can be used to obtain the cluster's total 
density distribution, which then allows one to determine the 
contribution due to SNs by subtracting the known density of 
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the ICM. Considering the equation of state for a partially de- 
generate neutrino gas, the resulting SN density p Y is then used 
to infer the associated radial velocity dispersion cr v needed 
for equilibrium. A detailed description of the actual calcula- 
tion can be found in Sec. 3 of Ref. Il83l . To check whether 
the results for p v and <x v are compatible with each other, one 
can exploit the TG phase-space constraint ||92ll . Assuming a 
Maxwellian velocity distribution, the maximally allowed den- 
sity Pv.max for a given value of <x v reads 



2 (2tt) 3 / 2 /z 3 



where the number of allowed helicity states is assumed as 
g v = 2 If83l IT331 and m v = lleV. For the "simple" MOND 
interpolating function which is defined as 



p(x) = 



1 + X 



(24) 



it has been found that the calculated SN phase-space density 
of all considered systems reaches the TG limit in the central 
part (r < 20kpc for A2390) [83), meaning that the SNs ac- 
quire their densest possible configuration in that region. If the 
equilibrium assumption is valid, this result further implies that 
a small portion of the dynamical mass must be covered by the 
brightest cluster galaxy. As for the cD galaxy of A2390 and its 



contribution in this context, we refer the reader to App. A 3 



In principle, we could directly adopt the SN density of 
A2390 calculated in Ref. [83] for our simple cluster model if 
we specified a TeVeS free function p which corresponds to the 
choice Eq. ( [24] > in MOND. For numerical reasons discussed in 
Ref. If59l (hereafter Paper I) and to maximize possible MON- 
Dian effects, however, we assume a TeVeS free function of the 
following form: 



(25) 



where y is defined according to Eq. Q. Apart from its sim- 
plicity, Eq. ( p5| > is close to Bekenstein's original choice of 
the free function (Paper I), and thus allows one to derive the 
TeVeS lens properties in a fully analytic way for certain con- 
figurations like, for example, spherically symmetric lens mod- 
els l56l . In the intermediate and low acceleration regime, 
which is typically realized in galaxy clusters, the MONDian 
counterpart of Eq. (|25j) can be expressed as |45 1 



p(x) 



Vl +4x- 1 
VI + 4x + 1 ' 



(26) 



which is known to yield a less favorable description for the 
rotation curves of spiral galaxies than Eq. tfZfy as it enhances 
gravity too efficiently [ 134|. Inserting the above into Eq. ( |22] i, 
we have repeated the analysis of Ref. |83| for A2390, and 
calculated the equilibrium SN density distribution suitable for 
our cluster model in TeVeS. The resulting density profile is 
shown as a dashed line in Fig. [2] Note that the apparent wavi- 
ness is not a numerical artifact, but rather emerges from using 
the data of Ref. Ml 081 in Eq. pi) . As the free function Eq. 



p6| ) enhances gravity more efficiently than Eq. ( |24| i, the SN 
density is notably decreased (cf. Fig. 2 of Ref. 1831 ). with the 
effect becoming stronger for larger radii. In the center, how- 
ever, there is basically no change, indicating that the previous 
constraints due to the TG limit remain the same. To simplify 
the input into a numerical solver, the obtained SN density can 
be well fit by a profile of the following form: 



P(r) = 



PO 



(1 + (r/r ) r ) 



(27) 



(23) where p ~ 5.5 x lO 7 M kpc 



-3 



14kpc and y ~ 8.2. For 



comparison to the numerical result (dashed line), the analytic 
fit (solid line) is also illustrated in Fig. [2] 

Note that the actual choice of the free function, which fixes 
the equilibrium distribution of SNs, will have no significant 
impact on the results for quasiequilibrium configurations pre- 
sented in Sec. [V] While the main cluster potential will almost 
be the same - it is exactly the same in case of spherical sym- 
metry - for different p, only the effects of substructure, e.g. 
the contribution of individual galaxies in A2390, should be 
affected by the particular form of the function p. Therefore, 
our decision to use Eq. ( |25] l will result in optimistic estimates 
of effects intrinsic to the framework of TeVeS. Since we are 
interested in the regime of strong lensing, however, we expect 
these differences to be rather mild. 



V. QUASIEQUILIBRIUM CONFIGURATIONS 

As a first approach, we shall investigate the strong lensing 
properties of quasiequilibrium configurations based upon sev- 



eral variations of the cluster model presented in Sec. IV Al- 
though these models do not provide a realistic description of 
the cluster's core region, their study will be extremely useful 
to explore intrinsic TeVeS effects and to see whether TeVeS 
offers alternative mechanisms - different from those in GR - 
which can produce straight images. For the sake of clarity, 
we discuss details on the used numerical tools and the basic 
simulation setup in App. [Bj 



A. Analysis of the TeVeS lensing maps 

Considering the previously introduced equilibrium model 
of the cluster, we are still left with substantial freedom regard- 
ing the galaxies' line-of-sight positions which are not con- 
strained by observations and may vary over the cluster's extent 
which we define by the model's cutoff radius R = IMpc in- 



troduced in App. A 1 Also, to account for nonsphericity of 
the cluster, we shall allow an additional ellipticity for the cen- 
tral density distribution (x-ray and SNs) which is solely mod- 
eled within the observed plane. Together with a respective PA, 
this gives a total of 7 free parameters for our simple model if 
we fix the galaxies' M/L ratios. As for the range of elliptici- 
ties, we choose a maximum corresponding to an axis ratio of 
bja ~ 0.7. Moving significantly beyond this threshold would 
cause a severe mismatch to x-ray observations I106l lll7ll . thus 
yielding a rather unrealistic cluster description. 
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FIG. 3. TeVeS lensing maps for a quasiequilibrium configuration: Shown are the resulting convergence k (left), shear modulus y (middle), 
and critical curves for a cluster model with e = 0.7, PA = 115°, mass model M\ and line-of-sight configuration B. The triangles indicate the 
observed position of the straight arc. 




6, [arcsec] A [arcsec] 6, [arcsec] 

FIG. 4. Same as Fig. [3] but assuming PA = 133° and mass model Mi. 



Modifications of the overall density profile along the line 
of sight have already been studied in Paper I: Varying the 
lens' extent between two extreme configurations, a disklike 
and a strongly "cigar-shaped" lens, can cause changes of 
up to 10 - 20% in the lensing maps as well as the critical 
curves. For realistic cluster models lying in a range between 
these extrema, however, this effect is expected to be less pro- 
nounced, typically accounting for deviations on the order of 
a few percent. Therefore, we shall ignore such modifications 
in this work. Also, since the straight arc's sources are located 
close in redshift space (the corresponding distances D s and 
D = D c {Dd s /D s only differ by roughly 3%), we restrict our- 
selves to a single source plane for our analysis. Unless oth- 
erwise stated, we will always work with a lens and source 
redshift of zi = 0.23 and z s — 1, respectively. 

For different plausible cluster configurations, we have 
found no quasiequilibrium model that is capable of produc- 
ing (nearly) straight images at the observed arc's position. As 
all of our results are qualitatively very similar, we only present 
a selection of simulation runs in the following. For example, 
Fig. [3]shows the calculated TeVeS lensing maps assuming an 
axis ratio bja « 0.7 corresponding to an ellipticity of e — 0.7, 
PA = 115°, mass model M\ and line-of-sight configuration B 
(see Table A similar case is illustrated in Fig. [4] assum- 



ing PA = 133° and mass model while keeping all other 
parameters the same. We note that the lensing properties in 
the arc's vicinity are almost entirely dominated by the clos- 
est galaxy. The structure of the critical curves (right panel) 
already reveals that such models will produce strongly bent 
images with respect to the galaxy 2592 at the position of in- 
terest. To elucidate this point and to demonstrate the problems 
of such configurations, we have constructed a luminosity dis- 
tribution, which roughly resembles the observed image mor- 
phology. This distribution has then been mapped back into the 
source plane, assuming the "second" cluster model presented 
in Fig. |4] The generated image and its associated source dis- 
tribution are shown in the left and right panels of Fig. [5] re- 
spectively. Our particular example exhibits several features 
indicating that the model is not compatible with the observed 
straight image. These features can be summarized as follows: 

(a) Around the area where the three source patches visible 
in the right panel of Fig. [5] appear to intersect (marked 
with a rectangle), the inferred source distribution becomes 
multivalued. This remains true even after taking into ac- 
count that the image is due to two distinct sources (see 
Sec. |III B| i, and thus the lens model turns out to be am- 
biguous and inconsistent. 
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FIG. 5. Left: The generated image contours (solid lines), resembling the observed luminosity distribution of the straight arc, and the critical 
curves (dashed lines) for an equilibrium cluster model with e = 0.7, PA = 133°, mass model M 2 and line-of-sight configuration B. Right: The 
resulting source distribution (solid lines) and lens caustics (dashed lines), where contours have been determined by averaging the calculated 
source points onto a regular grid. The open contour lines are due to a cutoff of the mapped image. In both panels, contours are in arbitrary 
units and chosen at equidistant levels. 




(b) Apart from the tangential caustic, i.e. the inner dashed line 
shown in the figure's right panel, the found source distri- 
bution also crosses the radial (outer) caustic, implying the 
existence of further images different from the straight arc. 
However, there is no evidence for such additional images 
as they are not observed in the system. 

(c) Assuming an average size of roughly 1" (~ lOkpc) for 
galaxies at z = 1, the source's constituents appear too 
big in angular size (up to 4"), yielding a rather unlikely 
scenario. This problem further deteriorates if one tries to 
avoid the issues related to (a) and (b) by lowering the total 
mass of the nearby (lens) galaxy 2592. 

Observations further indicate the presence of several faint 
elongated objects whose orientation is approximately the 
same as that of the arc, with a scatter of only a few degrees 
||99l . Together with the above, these arclets strongly support 
the requirement for a special lens composition rather than the 
necessity for unusual source properties, suggesting that the 
lens configurations considered here are inappropriate to ex- 
plain the straight image. 

This is the basic result of all simulated cluster models, 
which seems to be insensitive to the used mass model {M\ 
or Mj) or the actual line-of-sight alignment of galactic com- 
ponents. To quantify the effect of the latter, we compared the 
lensing maps of individual models for two extreme line-of- 
sight configurations, A and B. Adopting the parameters of 
the realization presented in Fig. [3] Fig. [6] displays the ob- 
tained relative difference between the corresponding conver- 
gence maps. As we can see, the deviation can reach values up 
to ~ 30% in regions of low (effective) surface density, but re- 
mains smaller (< 15%) in regions where k> \. A comparison 
of the corresponding critical curves and caustics of galaxies 
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FIG. 6. Relative difference of the convergence maps calculated for 
the two line-of-sight configurations A and B: The here presented re- 
sult assumes e = 0.7, PA = 133°, and mass model M\ . 

reveals that this line-of-sight effect typically affects their po- 
sition on the order of ~ 10%, which has no qualitative impact 
on our results. As for the dependence on the actually used 
mass models (galaxies), we will investigate the influence of 
varying M/L ratios in the next section. 

B. Variation of mass-to-light ratios 

So far, we have restricted our analysis to two sets of M/L 
ratios for the cluster's galaxies (see Table How robust are 
our results with respect to variations of these ratios? Here we 
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take a simplified approach to obtain reasonable estimates of 
the effect on the strong lensing properties, in particular, the 
critical curves. In what follows, we use a Hernquist profile 
H1321 with fixed core radius rn = 3kpc for the density dis- 
tribution of galaxies [corresponding to the limit e = in Eq. 



( |A2| >]. In the isolated case and for our choice of the free func- 
tion //, this allows one to express the lensing properties fully 
analytically (e.g., Paper I), and the deflection angle is given 
by " 



where 



VGMflo AGM% 



4GM£ 



(28) 



A(0 = 



(29) 



arsinh ^|l - (r H /£) 2 | € < r H 
arcsin ^1 - (r H /£) 2 g > r H 

Furthermore, we will assume that 



(a) the superposition principle remains valid; i.e. the lensing 
maps of isolated galaxies and the cluster background can 
just be added, which is rigorously true if the components 
are infinitely separated from each other and leads to an 
optimistic estimate otherwise, and that 

(b) the cluster background at each galaxy's position can be 
modeled as an external contribution with locally constant 
convergence kq and shear modulus jc- 

Choosing polar coordinates, the latter yields an effective clus- 
ter deflection potential of the following form: 

m ip) = K j-0 2 + ^6 2 cos(2( V - <p )), (30) 

where the external shear's principle axes system is defined by 
ipQ. Locally, the system's total shear modulus, relevant for the 
determination of critical curves and caustics, depends nonlin- 
early on the contribution due to the Hernquist lens and the 
cluster, 

r t 2 ot = (ri,n + ri,c) 2 + (yi,H + n,c) 2 - (31) 

Because of the shear's tensor property, the above is 
anisotropic, which directly affects the resulting position of 
critical curves given by Eq. (TT9J, 



(1-k) -y =0. 



(32) 



To obtain the mean effect due to the cluster background, we 
perform an average over <p and all possible orientations of the 
external shear field, which leads to 



t 2 2 

7-7h + Jc- 



(33) 



We use the simulation result for an equilibrium cluster model 
with e = 0.7, PA = 115° and no galaxies to estimate the pa- 
rameters of the background model. Around the arc's position, 
this roughly fixes kq ~ 0.29 and jc « 0.17. For this case, 
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FIG. 7. Predicted mean radius of tangential critical curves of a Hern- 
quist lens embedded into the external cluster field: The mean radius 
is plotted as a function of the projected enclosed mass within an aper- 
ture of 3" diameter (~ 1 lkpc). The vertical lines indicate the values 
of the galaxy 2592 for mass model Mi (dashed line) and (dotted 
line), respectively. 



Fig. [7] shows the resulting mean radius of the tangential crit- 
ical curve as a function of the enclosed galactic mass within 
an aperture of 3" diameter. While this should give a reason- 
able picture for the galaxy 2592, which resides close to the 
arc, the such estimated radii will be too large for the other 
galaxies. These are located in regions where the background 
has a weaker impact (k c and y c take lesser values), leading to 
an optimistic prediction of their mean critical-curve size. As- 
suming small variations in M/L, the figure suggests no qual- 
itative changes of our previous results. Even if we consider 
that M/L ratios may change up to a factor of 4 in the infrared, 
we find a maximum increase of the mean critical-curve ra- 
dius corresponding to a factor of approximately 2. At most, 
such an extreme scenario could come close to a merged-cusp 
model for the galaxies 2592 and 6666, but this configuration 
cannot explain the arc due to its inappropriate orientation and 
position of critical curves and caustics in the lens and source 
plane, respectively. We therefore infer that the basic result of 
Sec. V A does not depend on the particularly assumed M/L 

unless the background potential 



ratios of individual galaxies 
is substantially modified. 

We have also explored the influence of perturbations to the 
central cluster profile. For this reason, we have assumed a sec- 
ondary spherical clump made out of gas and SNs which fol- 
lows the same profile as the central distribution and accounts 
for 10 - 15% of the system's total mass. The clump's position 
has been chosen from a narrow range roughly centered on the 
detected substructure in the x-ray map 1 106 . 117] (see Table 
ttj. Again, we have found no qualitative difference compared 
to previous simulations. The calculated deviations in the lens- 
ing maps are on the order of a few percent, leaving a basically 
negligible impact on the critical curves and caustics. Similar 
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FIG. 8. Simulation results with proper treatment of the scalar field equation compared to those where the curl field h has been set to zero: The 
figure illustrates the relative deviation in the corresponding lensing maps, the convergence k, and the shear modulus y, assuming an equilibrium 
model with e = 0.7, PA = 133°, mass model M\ , and line-of-sight configuration B. The visible gridlike structure is a combined effect of Fourier 
fluctuations, interpolation, and the division by values close to zero. 



statements apply to an overall increase of the central density 
profile by 10 - 20%. 

Together with the results presented in Sec. 



VA we thus 



conclude that TeVeS quasiequilibrium configurations with 
lleV SNs are not capable of explaining the observed arc. In 
particular, we find no evidence for the formation of beak-to- 
beak or lips catastrophes [115] due to intrinsic TeVeS effects, 
which could give rise to straight images. Therefore - just as 
in GR - a suitable TeVeS lens model needs substantially more 
mass as well as a special density distribution in the cluster's 
core region. A general procedure on how to obtain such mod- 
els will be discussed in the next section. 



VI. NONEQUILIBRIUM CONFIGURATIONS 

In the following, we shall outline a general approach for 
modeling cluster lenses in TeVeS which allows one to use ex- 
isting GR lens models to estimate the needed TeVeS lens prop- 
erties. Adopting a bimodal lens model for the straight arc, we 
will present an example of such a lens and discuss implica- 
tions for the modified framework. 



A. Systematic approach to cluster lenses 

Taking a naive point of view, one might expect that strong 
lensing is subject to the strong acceleration regime, and there- 
fore it should be enough to consider the limit fi — > 1. In this 
case, all relevant equations would reduce to their GR coun- 
terparts, allowing a conventional lensing analysis. Previous 
calculations (Paper I) have shown that such an approximation 
is not justified. In particular, the scalar field can have a signif- 
icant impact on the second derivatives of the lensing potential. 
For instance, this can increase the radii of critical curves by up 



to a factor of 2, depending on the assumed mass distribution of 
the lens (cf. Figure 7 of Paper I). As we shall see below, how- 
ever, there is another way of simplifying the lensing problem 
in TeVeS. 

Let us return to the scalar field given by Eq. Q. Integrating 
once, we may recast this equation as 



An 



(V<& N + V xh), 



(34) 



where h is a regular vector field determined up to a gradient 
by the condition that the curl of the right-hand side of Eq. 
p4| ) must vanish. We note that the main difficulty associated 
with solving the scalar equation are the generally nonvanish- 
ing components of h. If for any reason h*0, Eq. p4| reduces 
to a relatively simple algebraic relation between the gradients 
of the scalar and the Newtonian potential, 



LlV(p ~ — V<l>Af, 

An 



(35) 



which can easily be inverted by numerical means to give V0, 
assuming that the Newtonian potential (or only its gradient) 
are known. Therefore, we want to address the question of 
how the field h is affecting the corresponding lensing maps 
in the strong lensing regime. We already expect h to be im- 
portant around local extrema of the Newtonian potential, but 
it is difficult to make any intuitive guesses about its quantita- 
tive impact in stronger gravity regions as well as on the final 
projected result. The most straightforward approach to this 
problem is a direct comparison of simulations treating the full 
scalar equation to those where h = 0. To this end, we have 
taken our previous quasiequilibrium models and fed them into 
a modified version of our solver, now assuming Eq. <j35j to de- 
termine the scalar gradient. Since our choice of fi is very close 
to that presented in Ref. l56l (e.g., Paper I), our code assumes 



|V0| = Vao|V<D A 



(36) 
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FIG. 9. Spherically averaged density (left panel) and corresponding projected surface density (right panel) profiles for the bimodal lens model 
in TeVeS: Shown are the results for the Newtonian dynamical mass (solid line) and the two TeVeS free interpolating functions corresponding to 
Eq. l |40[ > (dotted line) and Eq. i ]4 1 (dashed line), respectively. Note that the three-dimensional mass density profiles are entirely dominated by 
the contribution of SNs within a radius of a few hundred kpc. For r < 25 - 30kpc, the derived densities are well below the TG-limit-saturating 
equilibrium distribution of lleV SNs in Fig. [2] and they are also much broader. This already indicates that the main component's phase-space 
limit is not violated here. 



to calculate V</>. For instance, adopting an equilibrium model 
with e = 0.7, PA = 133°, mass model M\ and line-of-sight 
configuration B, the relative deviation of the lensing maps 
from the proper solution is presented in Fig. [8] While the 
convergence varies by 5 - 15%, differences in the shear map 
can be as high as ~ 50%. As expected, the largest deviations 
occur in regions where the Newtonian gradient approaches the 
null vector, which, for example, can be seen in the very core of 
the central elliptic profile for both the convergence and shear 
maps. 



B. Modeling the straight arc in TeVeS 



Clearly, the impact of the curl field h is not negligible in 
regions of low gravity. Concerning the domain of strong lens- 
ing, however, we find the following: Comparing the corre- 
sponding critical curves and caustics, the curl field turns out 
to be much less important. Interestingly, the obtained de- 
viation with respect to their position in the lens and source 
plane, respectively, is only about < 2 - 3%. Within a suffi- 
ciently large environment around these curves accounting for 
all strong lensing features, the accuracy of the approximated 
(h = 0) lensing maps is typically of the same order, meaning 
that the curl field negligibly contributes to the strong lensing 
properties of a given matter distribution. Our result appears 
to generally hold for strong cluster lenses and indicates that it 
is enough to consider Eq. ([35) in the context of TeVeS lens 
models. Therefore, if one specifies the line-of-sight extent of 
the total system as well as individual matter components, this 
offers a direct systematic way of modeling strong lenses in 
TeVeS. 



Based on the result of Sec. VIA one could, in principle, 



take an available GR lens fitting routine, modify it to include 
the TeVeS scalar field according to Eq. ( |35| l, and use it to ob- 
tain a lens model for the straight arc. It is obvious that such an 
approach will be computationally more demanding because 
the scalar's contribution has to be evaluated in three dimen- 
sions, and one also needs to invoke numerical integration to 
derive the desired projected quantities. In the following, how- 
ever, let us consider an alternative way to estimate the nec- 
essary deflection mass and its distribution in TeVeS. For this 
reason, we start from the bimodal GR model derived in Ref. 
Ill 171 which, in addition to the central matter clump, assumes 
a smaller subcomponent at approximately 45" (~ 166kpc) 
from the cluster center. The second clump is motivated by 
the existence of substructure in the cluster's x-ray map which 
is used to infer its position in the lens plane. Both clumps 
are chosen to follow a pseudoisothermal elliptic mass distri- 
bution (PIEMD) 1 135 1, but the subcomponent's profile is as- 
sumed to be spherically symmetric. Correcting for the here 
used cosmological background, the model gives an enclosed 
projected mass of M a ~ 1.2 x 10 14 M Q within a circular aper- 
ture of 38" (~ 140kpc) radius from the cluster center. As 
typical for strong lensing mass models, this estimate should 
lie within ~ 30% of the true value 11361 . 

Using the arguments presented in Sec. |VI A| it is obvious 
that there exists an analogous bimodal lens model in TeVeS. 
To obtain a spherically averaged density estimate in TeVeS, 
we ignore the secondary clump, which negligibly contributes 
to the enclosed mass within the given aperture, and also as- 
sume that the main component can be described by a spher- 
ically symmetric density profile. Thus, its three-dimensional 
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TABLE II. Fiducial parameters of the bimodal lens configuration pre- 
sented in Ref. 11171 : Here the subclump is offset by approximately 
45" (~ 166kpc) from the main component. 





b/a 


PA 


rc 


O"co 






o 


(") 


(kms- 1 ) 


Central main clump 


0.71 


49.2 


12 ± 5 


950 ± 100 


Subclump 


1 




7-12 


420 - 500 



matter distribution can be written as 



P(r) = Po 



r 2 c + r 2 ' 



(37) 



where po is the central density and rc the core radius. Alter- 
natively, Eq. ( f3"7] i may be written in terms of its asymptotic 
velocity dispersion cr^ associated with the density profile of a 
singular isothermal sphere: 



P(r) 



1 



2nG r 2 + r 2 



(38) 



The corresponding enclosed mass of this density distribution 
at radius r reads 



M(r) - 4-nr c po Ir - r c arctan 



rc 



(39) 



Since our choice of the free function allows us to make use of 
Eq. prj) , it is possible to express the enclosed mass in TeVeS, 
which effectively generates the same dynamical mass as Eq. 
(07), as 



Me ff (r) = M(r) + 



= M(r)- 



2G 



1- , 1 



4GM(r) 
a r 2 



4s 



(40) 



(l + VTT4^) 



where s = GM(r)/aor 2 . As previously noted, however, the 
choice of Eq. ( |25] l does not yield a good description of galaxy 
rotation curves. Adopting a TeVeS free function correspond- 
ing to Eq. fl24), a similar calculation leads to 



M eff (r) 



[M(r)Y 



M(r) + a r 2 /G 



M(r) 



1 + s 



(41) 



Setting r c * 13" (48kpc) 111 171 and requiring that the en- 
closed projected dynamical mass within 38" is still given by 
M a , the above expressions can be used to derive the under- 
lying density distributions which, together with the resulting 
surface density profiles, are illustrated in Fig. [9] The visi- 
ble density drop-off within r < 20kpc is a consequence of the 
assumed PIEMD and probably unphysical, but can easily be 
avoided by changing the central profile in favor of a peaked 
and finite core, fixing the enclosed mass around r = 140kpc 
(and thus keeping the lens properties needed for the arc). Of 
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FIG. 10. Estimated 1 leV SN density distribution (solid line) and 
corresponding TG limit (dashed line) for a subclump model with 
r c ~ 7" (26kpc) and cr^ « 500km s -1 : The origin is centered on 
the subcomponent, and the TG limit has been calculated according 
to Eq. ([23}, following the prescription of Ref. 1831 . Note that the 
slight "wiggly" feature of the dashed line is due to a nonuniform dis- 
persion a(r) which is computed in a self-consistent way 1 83 1. 



course, our results depend on the assumed line-of-sight extent 
specified by Eq. p7j ), but the derived surface densities should 



vary by only a few percent for different models (see Sec. V A i. 
We also note that the "modified" density profiles yield a finite 
mass; for both Eqs. ( |40) and ( |4T| ), the total mass is given by 
(taking the limit s — > 0) 



\6n l Gr c pl 



lim M eff (r) = 



(42) 



Although the profiles are not diverging, the relevant mass typi- 
cally extends to large radii. Therefore, the use of such profiles 
within the full TeVeS solver is very inefficient because very 
large box sizes would be necessary to perform the calcula- 
tions, which underlines the advantage of neglecting the curl 



field for strong lensing models (see Sec. VI A I. 

The resulting total lensing mass is entirely dominated by 
SNs within a radius of a few hundred kpc, which allows 
one to ignore the contribution of gas and stellar material to 
excellent approximation. We have checked that the three- 
dimensional density distributions in Fig. [9] basically rep- 
resenting the lens model's main component, are consistent 
with the TG limit estimated for hydrostatic equilibrium and 
a Maxwellian velocity distribution, following the approach of 
Ref. l83l . This is already indicated by the fact that the de- 
rived densities are much broader and well below the TG-limit- 
saturating 1 leV SN equilibrium distribution (shown in Fig. [2]) 
for r < 25 - 30kpc. At the arc's position (6 = 38"), the actual 
enclosed projected mass of the TeVeS lens models is given by 
6.1 x 10 13 M o or 8.0 x 10 13 M G , assuming Eq. gO) or tgl), re- 
spectively. Here the model's subcomponent deserves special 
attention: Naively treating the problem, the smaller clump's 
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the central density which could affect our estimate of the TG 
bound. To check this, we have modified the central SN den- 
sity profile of the subclump model presented in Fig 
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FIG. 11. Same as Fig. [10] but assuming r c as 10" (37kpc) and 
cr„ as 440km s -1 . 



presence acts as a perturbation to the total system's phase- 
space density, and thus it is trivially in accordance with the 
estimated TG limit since the main clump is. However, this ap- 
proach typically leads to overestimating the TG limit, consid- 
ering that the secondary clump should be regarded as a bound 
object by itself. Taking the view that A2390 has undergone 
recent merger activity, it seems reasonable to assume that the 
subcomponent has formed at a sufficiently earlier time, and 
therefore it should be subject to its own phase-space distri- 
bution. This suggests that one should examine the secondary 
clump separately. 

Considering the subclump as an isolated object entirely 
dominated by 1 leV SNs, we have repeated the above TG anal- 
ysis for the range of parameters listed in Ref. Ill 171 (see Ta- 
ble [n). To achieve a rather realistic TeVeS mass estimate, we 
have adopted a free interpolating function corresponding to 
Eq. ( [24] > for our calculations. The obtained SN density pro- 
file and the TG limit according to Eq. ( |2"3) are illustrated for 
two cases in Figs. 



10 and[TT 



Assuming rc ~ 7" (26kpc) 
and (Too ~ 500km s _1 , the subcomponent's density slightly 
exceeds the TG limit (up to 30%) within a range of approxi- 
mately 10 - 25kpc. Moving toward larger radii (r > 50kpc), 
the SN density consistently stays below this limit. For a less 
compact model with r c as 10" (37kpc) and <x M as 440km s , 
the TG bound is never exceeded. Generally, our results seem 
to rule out configurations where the subclump is modeled with 
small values of rc (< 8 - 9") whereas the bimodal TeVeS lens 
appears consistent with lleV SN HDM for larger choices of 
the core radius. Before drawing such a conclusion, however, 
we need to consider how strong the implication of the present 
analysis really is. 

First of all, we note that the lens model is based on the 
PIEMD model given by Eq. (37) . This clearly introduces a 
bias on our estimates; other assumptions about the compo- 
nents' density distributions might yield a different result. In 
particular, the PIEMD model leads to an unphysical drop of 



favor of a uniform core, but without changing its properties 
beyond r ~ 15kpc (the arc appears at r as 26 from the subcom- 
ponent's center). The resulting density profile and the corre- 



sponding TG limit are shown in Fig. 12 While the TG bound 
is still violated within ~ 10 - 25kpc, we see that the density 
limit is notably decreased in the center, almost matching the 
assumed SN distribution. Therefore, it is unlikely that shift- 
ing matter to the central region can help to avoid an excess of 
the TG bound. Next, our estimates assume that the subclump 
can be treated as an isolated object. Since the clump resides 
within the background field of the main component, this is not 
rigorously true. Using Eq. pT) , we find that the main com- 
ponent provides an external Newtonian gravitational field of 
around ao at the subclump's position. As for the subcompo- 
nent, this modifies the relation between gravitational field and 
underlying density distribution, and gives rise to an increase 
of the central SN density on the order of unity. Such a den- 
sity boost could push seemingly consistent subclump models 
with r c > 10" toward or even beyond the TG limit, but de- 
tailed statements about this issue are very sensitive to the ac- 
tual model parameters. 

Another point is related to the fact that our calculations 
rely on completely SN dominated lens components within 
~ lOOkpc. If placed at or close to the subclump's center, al- 
ready a relatively small, concentrated baryonic mass, e.g. a 
galaxy, on the order 10 8 - 1O 9 M could help to relax the den- 
sity constraint due to the TG limit 11371 . Whether such an 
approach can be reconciled with observations of this region, 
however, remains to be seen. Last but not least, we also need 
to check the viability of the current estimate of the TG limit 
which has been derived under simplified conditions. In what 
follows, we shall discuss in more detail how these simplifica- 
tions affect our analysis. 

As previously mentioned, the strong lensing domain in the 
center of A2390 is not in equilibrium and has a rather compli- 
cated nonsymmetric density distribution. This will obviously 
have an influence on the estimates for the TG limit. Consid- 
ering nonequilibrium configurations in general, the velocity 
dispersion cr is expected to increase for a given matter distri- 
bution when moving away from equilibrium. Since the value 
of cr increases in this case, one would also obtain a higher 
density limit for SNs according to Eq. ( |23) . Taking the ad- 
ditional asymmetry into account, however, the situation be- 
comes less clear. Depending on the system's actual proper- 
ties, the TG limit could increase or decrease, and there seems 
to be no universal rule allowing one to make solid statements 
for anisotropic systems. Finally, one should also adhere to de- 
viations from a Maxwellian velocity distribution. This issue 
has been addressed in Ref. 111381 . There it has been shown 
that the actual physical density limit becomes larger than the 
previous estimates of the TG limit which can be exceeded by 
up to a factor of 2. Again, this would imply that SNs could 
account for more mass in the cluster. Combining the three as- 
pects from above, it seems reasonable to assume that the true 
density limit will be on average higher than our previous es- 
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FIG. 12. Same as Fig. [TO] but now assuming a uniform core of the 
SN distribution. Note that the radius at which the density becomes 
constant (r as 15kpc) is fixed by requiring a continuous distribution 
ofSNs. 



timates, meaning that density models with SNs become more 
flexible. Note, however, that such an argument generally does 
not replace the need for a rigorous treatment of particular sys- 
tems. 

Given the accuracy of our present analysis and account- 
ing for all of the above, we conclude that the bimodal TeVeS 
lens model for the arc is in accordance with the assumption of 
lleV SNs. Nevertheless, it seems intriguing that the needed 
amount and distribution of 1 leV SNs lies so close to what they 
can maximally contribute to the system. It should be obvious 
that all of our statements depend on the assumed lens configu- 
ration and are valid only for the bimodal model we have con- 
sidered here. In particular, the bimodal lens model ignores the 
contribution of galaxies. These can have a significant impact 



on the lensing maps (see Sec. V A i, which is especially true 



for the galaxy 2592 adjacent to the straight arc. A more real- 
istic approach including all visible components would be use- 
ful to further constrain the properties of additional substruc- 
ture and check whether such configurations remain consistent 
with respect to the TG limit. While our analysis is concerned 
only with the straight arc, the cluster A2390 actually exhibits 
a number of lensing features which should all be taken into ac- 
count for a complete cluster model. Extending the investiga- 
tion also to other massive galaxy clusters, future work should 
address such complex lens models and their implications for 
TeVeS or related theories and 1 leV SNs; a systematic way for 
approaching this problem has been outlined above. 



VII. CONCLUSIONS 

In this work, we have suggested the use of strong gravi- 
tational lensing by galaxy clusters as a test of the combined 
framework of TeVeS and massive SNs. Originally motivated 



by theoretical and recently also experimental particle physics 
E3 H39l [HD, the idea of SNs with a mass around 1 leV has 
gained further interest as it provides a possible remedy for the 
problems of TeVeS and related theories ranging from large 
cosmological scales down to galaxy clusters. Unlike conven- 
tional CDM, such a fermionic HDM component is subject to 
strong phase-space constraints imposed by the TG limit. This 
allows one to check cluster lens models inferred within the 
above framework (or related ones) for consistency. 

As an example, we have studied the cluster lens A2390 with 
its notorious straight arc. Because of its elongation and orien- 
tation, the straight image appears to be quite unusual and indi- 
cates the need for a rather special lens configuration. Adopt- 
ing the approximation for weak fields and quasistatic systems, 
one of the main problems associated with the lensing analysis 
is the nonlinear relation between the TeVeS metric potential 
and the underlying matter density distribution. This nonlin- 
earity prevents one from working with projected quantities 
and requires one to perform all calculations in three dimen- 
sions. In addition, one is left with a nontrivial, Poisson-type 
partial differential equation for the TeVeS scalar field. 

To make some progress, we have considered a class of 
cluster models, based on the assumption of hydrostatic equi- 
librium, and investigated their lensing properties. This has 
been achieved by employing a MPI parallel solver for the 
TeVeS scalar field equation and simulating the correspond- 
ing lensing maps on the HUYGENS supercomputer which is 
located in Amsterdam. Our results imply that such quasiequi- 
librium configurations are not capable of explaining the ob- 
served straight arc. In particular, we have found no evidence 
for the formation of beak-to-beak or lips catastrophes Ml 1511 
due to intrinsic TeVeS effects, which could give rise to straight 
images. Line-of-sight effects and the impact of perturbations 
are typically small, changing the quantities of interest only on 
the order of a few percent. Similar to the situation in GR, a 
suitable TeVeS lens model therefore needs substantially more 
mass as well as a special density distribution in the cluster's 
core region. 

Based on the above results, we have further outlined a gen- 
eral and systematic approach to cluster lenses which signifi- 
cantly reduces the problem's complexity by avoiding the need 
of solving the TeVeS scalar field equation. Combined with 
conventional lensing tools, this opens a new window to strong 
gravitational lensing in TeVeS-like modified gravity theories. 
As a first application, we have explored the TeVeS analog of 
the bimodal lens configuration discussed in Ref. 11 171 . For 
this model, we have derived the SN distribution necessary to 
produce the desired image, using a simplified approach. The 
obtained SN density profile has then been compared to the 
maximally allowed contribution set by the TG phase-space 
constraint. To this end, we have estimated the maximal den- 
sity due to the TG limit following the prescription of Ref. Il83l 
and found a slight excess of this limit for the model's sec- 
ondary component if its core radius is small (jc % 8" - 9"). 
For less compact models, however, the TG bound is not vio- 
lated. Given the accuracy of our current analysis, we therefore 
conclude that the bimodal TeVeS lens model appears consis- 
tent with the hypothesis of 1 leV SNs. 
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Note that the bimodal lens model ignores the contribution 
of galaxies. As has been shown, these can have a significant 
impact on the lensing maps. A more realistic approach, in- 
cluding all visible components and other lensing constraints, 
should be taken into account to obtain better bounds on the 
required SN distribution and to check whether such configu- 
rations remain consistent with respect to the TG limit. Fu- 
ture work should address more accurate ways of estimating 
the TG limit in this context, and we suggest extending the 
investigation to other massive galaxy clusters which indicate 
the need for dark substructure. Unless one considers different 
solutions to the missing mass problem inherent to this par- 
ticular kind of modifications (see, e.g., Ref. [18]), the basic 
approach presented here should apply to any class of tensor- 
vector or tensor-vector-scalar theory which recovers the dy- 
namics of MOND in the nonrelativistic limit. Lensing by 
galaxy clusters could therefore provide an interesting discrim- 
inator between CDM and such modified gravity scenarios sup- 
plemented by SNs. In addition to the above, we note that 
next-generation neutrino experiments I14U41431 will further 
constrain the plausibility of lleV SNs. Even if they remain 
viable candidates, it still needs to be seen whether such SNs 
do actually cluster in the desired way fl83j. 

Finally, we advert to the fact that our analysis neglects pos- 
sible contributions due to perturbations of the TeVeS vector 
field Such contributions are known to be crucial for the 
formation of large-scale structure iBOllBTl . where they provide 
the key to enhanced growth while perturbations of the scalar 
<f> only play a subordinate role. As already pointed out in the 
literature ll93ll . this typically affects scales > 0.1 - IMpc and 
could be important for galaxy clusters. Owing to the more 
sophisticated structure of the field equations, however, even 
a rough magnitude of the vector's impact on these scales has 
not been estimated yet. Thus our work emphasizes the need 
for a quantitative description of these vector instabilities on 
small to intermediate scales, i.e. ~ 0.01 - IMpc. We also note 
that the result of such an analysis could strongly depend on 
the particularly assumed theory. 

Despite the previously mentioned limitations of the present 
work, our numerical simulations are probably by far the most 
detailed in the context of TeVeS and certainly provide the first 
extensive study of strong lensing features within this modi- 
fied gravity framework. Applications of our grid-based lens- 
ing code (e.g. with respect to offsets between visible mat- 
ter and weak or strong lensing features ||9T1 1144 1) hold the 
promise of very constraining limits on TeVeS(-like) theory 
combined with HDM and other unified recipes for the dynam- 
ics of MOND and DM flTS] [T45lfl47l . 
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Appendix A: Modeling the baryonic content of A2390 

1. X-ray gas and central mass distribution 

To derive a reasonable model for the gas distribution in 
A2390, we use the results given in Ref. H106H . The intrin- 
sic electron density derived from CHANDRA observations 
(shown in Fig. 10 of Ref. [106|) can be well described by 
a spherical profile of the following form: 



n e {r) 



(l+Wr ) 2 ) 



-1/2 > 



(Al) 



where «o = 0.1cm -3 and ro = lOkpc. Assuming a mean 
molecular weight of w — 0.6 and an additional factor of 1.2 
to account for the global effect of the cluster's stellar compo- 
nents, we thus obtain an expression for the effective central 
density profile with a central density of po = 1.8 x 1O 6 M 
kpc~ 3 . Since the volume integral of Eq. (JATJ diverges, we 
smoothly cut the profile at radius R within a range of 200kpc. 
The cutoff scale is set to R = IMpc which corresponds to 
0.7r5oo 11481 as given in Ref. |108|. This yields a total inte- 
grated mass of M ~ 1 .3 x 1O 14 M and a surface density profile 
which is in good agreement with a 10 - 20% gas fraction of 
the enclosed projected lensing mass estimated in the frame- 
work of GR M1101I1121 . The density distribution specified by 
Eq. < | A 1 1 > is illustrated in Fig. |2](dotted line). 

Although our choice for the density profile is less accurate 
and results in a slightly smaller mass than typical £1 models 
M106II1 171 or more flexible ones H108L it will be sufficient for 
our analysis. As is shown in Sec. |IV B| the relevant lensing 
mass is mostly dominated by the contribution of SNs. Thus 
the strong lensing results, which we are primarily interested in 
here, will be relatively insensitive to the actual assumption of 
the central baryonic distribution. Adopting the more realistic 
density models above in a few selected simulation runs, other- 
wise identical to those presented in Sec. [V] we find only small 
differences on the order of a few percent in the corresponding 
results and confirm our argument. This is also indicated by 
comparing the enclosed projected dynamical mass profiles of 
our cluster model (gas + SNs) to the Navarro-Frenk- White 
(NFW) profile [149) estimated in Ref. ED (see Fig. [13 1. 
Although the TeVeS model underestimates the mass, the dis- 
crepancy from the NFW model is only about 10% at the arc's 
position (9 m 38")- In addition, the figure shows the weak 
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FIG. 13. Enclosed projected (Newtonian) dynamical mass profiles 
for our TeVeS equilibrium model (gas + SNs; solid line) and an NFW 
model (dashed line). The lower mass of the TeVeS model is mostly 
caused by the approximate description of the gas density given in 
Eq. ((ATJi; triangles indicate the estimates from weak lensing ob- 
servations. At the arc's position (9 « 38"), the relative difference 
between the models is about 10%. 



lensing results obtained from the Canada-France-Hawaii Tele- 
scope (CFHT) for a photometric redshift distribution based 
on the CFHT Legacy Survey data 11 121 11501 . The relative 
good agreement between dynamical and weak lensing mass 
estimates further implies that structure along the line of sight 
plays no significant role and does not affect our analysis. All 
presented quantities have been corrected for the cosmological 
model specified in Sec. 



II B 



Note, however, that a rather ac- 
curate description of the gas density as well as its temperature 
profile is important to estimate the neutrino content necessary 
for hydrostatic equilibrium in TeVeS 



2. Galaxy morphology and masses 

Since a rather detailed model of the cluster might be impor- 
tant in TeVeS, we also need to take the contribution of indi- 
vidual galaxies into account. For simplicity, we only consider 
the most massive member galaxies in the immediate vicin- 
ity of the straight arc's observed position; galaxies which are 
located farther away are unlikely to affect the TeVeS lensing 
maps in this region, which is confirmed by our results pre- 
sented in Sec. [V] Although A2390 exhibits a rich class of 
galaxy morphologies, with many galaxies showing elliptical 
or lenticular shapes, the impact of individual morphologies on 
the arc's environment can safely be neglected due to the galax- 
ies' sufficiently large distances. While this is not necessarily 
true for the galaxy 2592 (see Fig. [TJ which resides directly 
adjacent to the arc, a spherical density model provides a good 
description, which is indicated by the rather mild ellipticity 
seen in the optical HST image. As can be seen in Sec. [V] this 



approximation does not affect the basic results of our analysis 
- at least in the case of quasiequilibrium configurations. 

Furthermore, we assume that all considered galaxies can be 
modeled by a matter distribution of the form 



P(r) = 



2n(r + e)(r + r H ) 



3 ' 



(A2) 



where p(0) = M/(2ner 2 H ) is the central matter density, and 
the profile's core radius is universally set to r# = 3kpc. The 
length scale e corresponds to a smoothing parameter becom- 
ing necessary due to the limited resolution of our simulations 
and is specified in App. B 2 For e = 0, Eq. ( |A2[ ) reduces to 
the well-known Hernquist profile 013211 which closely approx- 



imates the de Vaucouleurs R 1 ^ law for elliptical galaxies. 

To infer the masses of individual galaxies, needed for our 
strong lensing analysis, we consider the data of the spectro- 
photometric catalog compiled in Ref. 11231 . which lists mag- 
nitudes for 48 galaxies inside the cluster A2390. All magni- 
tudes are given in the Gunn r band 11511 . and a simple formula 
1 152 1 to convert the R Johnson magnitude and the B—V color 
index to the Gunn r band can be found in the literature [ 153 1. 
Accordingly, we have computed r B , the Gunn r magnitude of 
the sun, adopting R Q = 4.42 fl54l and (B - V) Q = 0.64 lfT55l . 
We have found r Q = 4.95 which is rather close to the r value 
inferred from SDSS, the corresponding band being quite sim- 
ilar to the Gunn r band. Our result for r Q has then been used 
to evaluate the absolute luminosities of the galaxies given in 

Ref. nm. 

Next, we need a realistic mass-to-light ratio (M/L) in order 
to determine the galaxy masses. To this end, we have fol- 
lowed a twofold approach: First, we have adopted a constant 
M/L derived by combining the relation between M/L and the 
g — r color index presented in Ref. [ 156 1 with the g — r colors 
for massive ellipticals in the red sequence of the SDSS given 
in Ref. II 1 571 . The corresponding masses are labeled as M\. 
Second, we have also considered M/L as a function of M in 
agreement with the results for the galaxies of A2390 discussed 
in Ref. H123H . For this, a dynamical mass estimate based on 
measured velocity dispersions was used. As elliptical galax- 
ies are mostly subject to the strong gravity regime within their 
half-light radius, however, estimates in both MOND/TeVeS 
and Newtonian dynamics should be roughly the same. This 
second mass estimate, denoted as M%, is probably more reli- 
able since it involves fewer assumptions. The properties of the 
such obtained galaxy models are listed in Table [I] 



3. Role of the central cD galaxy 

Assuming an equilibrium model for A2390, it has been 
found that lleV SNs reach their densest possible configura- 
tion for r < 20kpc (see Sec. |IV B| i. Since Eq. ( p7| ) takes 
the TG bound into account, our cluster model misses some 
mass in the central part and does not correspond to a gen- 
uine equilibrium situation. A way of compensating for this 
is to consider an additional contribution due to the central cD 
galaxy. Following the lines of Ref. [83], one can estimate a 
total galaxy mass of approximately M - 1.8 x 1O 12 M . As 
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the central region of A2390 is neither spherically symmetric 
nor in equilibrium 1106111081 . it is important to note that such 
an approach has no real physical meaning, but rather offers a 
convenient way to tweak our cluster model. 

What does the above mean for our lensing analysis? Mod- 
eling the cD galaxy as a point mass, a straightforward calcu- 
lation shows that its impact on the TeVeS lensing maps can 
be safely neglected. At the position of the straight arc (38" 
or 140kpc from the cluster center), the additional matter gives 
rise to changes of 1 - 2%. Moving to smaller radii, the devi- 
ation grows, but we are not interested in this region anyway. 
Thus we consider the cluster model presented in Sec 
sufficient for our investigation. 



IV 



as 




Appendix B: Numerical tools and setup 

1. Solving the scalar field equation 

Having set the framework of gravitational lensing and cos- 
mology in Sec. [II] we may proceed with calculating the de- 
sired TeVeS lensing maps. The main problem associated with 
this task is to solve the scalar field equation specified in Eq. 
<j3j which can be rewritten as 

A0=p, (Bl) 

where the effective density p(p, d\<p, didjCp) is 

P ^ k J P ~ lk j% (09i<W)0W)) ' ( B2 ) 

and indices run from 1 to 3. Equation ( |B1| | corresponds to 
a nonlinear second order elliptic boundary value problem and 
can be tackled numerically. A Fourier-based solver employing 
an iterative relaxation scheme and an equidistant grid has been 
presented in Paper I where the basic algorithm and involved 
approximations are extensively discussed. 

One of the numerical challenges of our analysis is that we 
need to resolve galactic scales in a cluster-wide box, which 
requires a relatively large number of grid points. Since all 
calculations have to be performed in three dimensions, this 
clearly exceeds the capacity of a single-processor machine, in 
terms of both needed time and memory, and therefore calls 
for a more powerful computer architecture. For this rea- 
son, we have implemented a parallel version of the original 
solver using the Message Passing Interface (MPI) standard. 
The parallelization as well as all calculations presented in this 
work have been carried out on the HUYGENS supercomputer 
at SARA in Amsterdam within the HPC-EUROPA Transna- 
tional Access Programme. The HUYGENS system consists 
of 104 nodes, with 16 dual core processors (IBM Power6, 4.7 
GHz) as well as either 128 GBytes or 256 GBytes of mem- 
ory per node, thus providing an excellent environment for our 
needs. 

The parallel solver has been tested with analytic TeVeS 
models such as the Hernquist lens (see, e.g., Ref. |56| or Pa- 
per I), and has also been compared to previous calculations for 
the "bullet cluster" (Paper I), yielding exactly the same results 
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FIG. 14. Predicted critical curves for an isolated galaxy given by 
Eq. I |A2} : Assuming an aperture mass of 3.5 x 10"M G within a 3" 
(~ 1 lkpc) diameter as well as a lens and source redshift of zi = 0.23 
and z s = 1 , respectively, we present results for both a high resolution 
(~ 0.05", solid line) and a low resolution setting (~ 1.2", dashed 
line) with subsequent interpolation. 

- up to machine accuracy - as the serial version for identical 
input parameters. Considering the numerical setup for A2390, 
we choose a physical box size of V = cfi — (4Mpc) 3 in order 
to meet the requirements of the point lens approximation at 
the grid's boundaries (for details see Paper I). Performing a 
variety of test runs, we have found that the solver's conver- 
gence property quickly deteriorates if we increase the number 
of grid points per dimension N, meaning that the code takes 
many iteration steps or even fails to converge |158|. Typi- 
cally, this problem already occurs at N = 512 and manifests 
itself through extreme fine-tuning of the constant relaxation 
parameter a> which determines the scalar's mixing of at each 
iteration step, 

A0 w =p (n) , (n+1) = ojft n) + (1 - a))(f> (n \ (B3) 

Depending on the particularly used density model of the clus- 
ter, acceptable values for a> vary within a range of 0.7 - 0.9, 
but allow them to be easily identified just after a few iterations. 
Compared to the analysis of Paper I, we thus obtain no univer- 
sal value for the relaxation parameter. Similarly, we also note 
that the solver's behavior becomes more sensitive with respect 
to the scalar's initial guess. This is expected because the ef- 
fective deviation from the desired solution increases with N, 
and can usually be accounted for by slightly modifying the 
original point mass ansatz of Paper I to achieve a finite core, 

(O) (r) oc log(r + r c ), (B4) 

where r c is on the order of a few d/N. While more elaborated 
guesses are also possible, they typically do not yield a much 
better performance. 
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2. Numerical setup for A2390 

In all simulation runs, we set the number of grid points per 
dimension to N = 896. This yields a resolution of approx- 
imately 1.2" (~ 4.5kpc) for our choice of d — 4Mpc. To 
improve the numerical stability of our Fourier solver, we fur- 
ther require all density components to be centered within their 
respective subcube, which can lead to a maximal deviation of 
0.6" from the positions listed in Table [I] In addition, we as- 
sume a smoothing parameter e = lkpc for the galaxy profile 
given by Eq. ( |A2[ ). Once the desired fields and derivatives are 
calculated, we use a cubic spline to interpolate our results and 
determine the relevant lensing quantities. For the given speci- 
fications, individual simulation runs typically require 30 - 50 
iteration steps to converge, and can last up to 24 hours using 
32 processors. 

The interpolation approach is justified because the exact re- 
sult is expected to be relatively smooth. To support this argu- 



ment, we performed a small numerical experiment: Assum- 
ing an aperture mass of 3.5 x 1O U M within a 3" (~ llkpc) 
diameter and the parameters from above, we compared the 
predicted critical curves of an isolated galaxy given by Eq. 
( |A2| > for low resolution (~ 1 .2") with subsequent interpolation 
to those calculated for a higher resolution setting (~ 0.05"). 
Choosing a lens and source redshift of zi = 0.23 and z s = 1, 



respectively, the results are shown in Fig. 14 While the radial 
critical curve is not very well recovered, the radius of the tan- 
gential critical curve, which is relevant for our considerations 
on the straight arc II 1591 . is only underestimated by roughly 
10% on average. Considering the full cluster model of A2390, 
however, galaxies are not isolated, but reside within the clus- 
ter's background field, which leads to a boost of their corre- 
sponding Einstein radii. Therefore, we expect the accuracy 
of the calculated lensing properties, including critical curves 
and caustics, to be significantly improved and sufficient for 
our analysis in this case. 
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